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ABSTRACT 

The  ability  to  measure  and  predict  atmospheric 
turbulence  affecting  laser  baam  propagation  is  a  major 
concern  when  considering  military  applications.  Suoh  a 
method  using  a  telesoope,  high  resolution  television  camera, 
video  tape  recorder,  digital  storage  oscilloscope,  and 
calculator  system  has  been  devised,  tested  and  utilized.  A 
laser  beam  signal  is  recorded  on  video  tape  for  further 
processing.  This  signal  is  displayed,  stored  and  digitized 
using  a  Tektronix  453  digital  storage  oscilloscope.  The 
digitized  signal  is  ssnt  to  a  Hewlett-Packard  9825  computing 
system  for  Fourier  transform  analysis  and  determination  of 
the  refractive  index  structure  constant,  ^n2.  Several 
trials  were  conducts!  using  He-Ne  and  Ga-As  lasers.  The 
results  demonstrate!  good  correlation  with  theoretical 
predictions  as  well  as  previously  analyze!  lata. 
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I.   INTRODUCTION 

A.       BACKGROUND 

The  increasing  use  of  lasers  and  laser  technology  for 
military  applications  has  brougit  about  a  need  for  analysis 
of  the  laser  beam  in  its  environment,  the  turbulent 
atmosphere.  A  projast  at  Naval  Postgraduate  School  dealing 
with  this  subject  is  continuing  and  is  the  main  topic  of 
This   thesis. 

The  patterns  produced  by  lasers  on  targets  have  inherent 
problems  that  include  broadening,  beam  wander  and  intensity 
fluctuations  brought  about  by  turbulence  in  the  atmosphere. 
These  effects  of  atmospheric  turbulence  on  laser  propagation 
have  been  well  deternined  [  Ref .  1].  In  terms  of  the  Fried 
model,    [Ref.   2],  C  n2,      the      refractive      index      structure 

constant,  has  been  determinel  to  adequately  express  the 
Modulation  Transfer  Function  (MTF) ,  or  the  Mutual  Coherence 
Function    (MCF)     for    the    atmosphere. 

A  system  has  been  developed  that  provides  a  measurement 
of  C  n  2  for  atmospheric  turbulence  along  the  optical  path 
through    which    a    laser    is      propagated.         The    system    employs   a 


vidicon  and  telescope  as  the  detector  and  a  distant  laser  as 
the  source.  This  system  duplicates  the  slit  scanning  system 
presently  in  use  at  Naval  Postgraduate  School  [Ref.  2]. 
Tests  using  the  vidicon  equipment  have  been  previously 
completed  with  values  for  Cn2  on  the  order  of  10~15  m~2/3 
being  obtained  [Ref.  3]..  Measurements  have  been  made  using 
the  same  experimental  set-up,  sxcept  that  the  Tektronix  468 
Digital  Storage  Oscilloscope  are  used  for  digitization 
instead  of  the  Quantise  DS-30  Digital  Video  Processor.  These 
measurements  have  demonstrated  values  of  Z n2  of  comparable 
accuracy. 

B.        PROBLEM 

By  using  the  digital  storage  capability  of  the  Tektronix 
468,  the  previously  recorded  laser  signal  is  used  as  an 
input  and  is  evaluated  by  modifying  the  program  developed  by 
Crager  [Ref.  4].  a  brief  overview  of  the  procedure  is 
described  below.  Detailed  explanticns  of  the  experiment  and 
computer  program  are  contained  in  Chapters  III  and  IV, 
respectively. 

The  approach  taken  is  basically  the  same  as  that  of 
Crager,  but  because  of  the  unavailability  of  a  disk  ROM  it 
is   necessary   to   store   both  raw   and   processed   data   on 
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magnetic  taps.  Tha  data  a:a  loaded  and  recorded  as 
necessary  by  the  HP  33  25.  The  basic  assumption  permitting 
analysis  is  that  the  horizontal  TV  scan  line  through  the 
laser  spot  is  considered  to  accurately  mirror  a  point  spread 
function  of  the  inage.  Final  analysis  using  the  468 
oscilloscope  has  shown  this  to  ba  a  valid  hypothesis.  The 
digitized  data  from  the  468  3  scilioscope  agree  with  the 
previously  measured  lata  which  iss  many  pixels  of  television 
data  digitized  by  tie  DS-30  to  express  the  point  spread 
function. 

The  sequence  of  analysis  is  that  a  video  tape  recording 
is  made  of  the  TV  image  of  the  propagated  laser  beam.  The 
output  of  the  video  recorder  is  sent  to  che  Tektronix  468, 
where  the  derived  TV  scan  line  is  digitized,  averaged,  and 
stored.  The  HP  9825  records  the  digitized  data  and  produces 
a  line  spread  function  (LSF) ,  by  integrating  the  point 
spread  function, 

y  (max) 
L3F(x)  =   /pSF(r)  dy  (1.1) 

y  (min) 

where 
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E    =    yJ(*2+Y2) 


1 


and      computes      the      Fourier      transform      of      the      LSF.  The 

diffraction  limited  Fsurisr  transform  of  the  optics  is  now 
computed  if  the  MTF  of  the  optics  has  not  been  previously 
measured  experimentally.  Next,  the  program  finds  the  MTF  of 
the  atmosphere  by  dividing  tha  Fourier  transform  of  tha  LSF 
by  the  Fourier  transform  of  tha  optics.  Finally,  by  curve 
fitting,    the  program   computes    a    single   value    for   Cn2. 

The  program  now  predicts  tha  size  of  a  laser  spot  on  a 
target  using  the  calculated  valaa  of  Cn2.  This  value  of  Cn2 
is  used  to  calculate  an  MTF  of  the  atmosphere  which  is  then 
multiplied  by  the  Fourier  transform  of  a  source,  and  the 
Fourier     transform        Df      the      optics.  The         program      then 

calculates  the  inverse  Fourisr  transform  of  the  products  of 
the  above  and  uses  the  Abel  transform  to  give  the  angular 
point    spread  intensity    distriDution.  From    this   data,      the 

fraction  of  energy  as  a  function  of  the  total  energy  within 
a   given   radius    R   is    calculated. 
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II.   THEORETICAL  Z 3NSIDERATI DNS 

The  theory  of  laser  beam  propagation  through  a  turbulent 
medium  has  been  explained  by  Crittenden,  and  others,  and  is 
re-emphasized  here  for  continuity  purposes  [Ref.  2].  Since 
the  effects  of  turbulencs  on  laser  weapons  is  of  major 
concern,  measurement  and  prediction  capability  for  these 
effects  on  laser  beams  is  entirely  relevant. 

Due  to  the  existence  of  reciprocity,  these  ideas  apply 
to  either  laser  designators  or  imagers  [Ref.  5].  In  the 
study  of  the  theoretical  modeL  by  Fried,  the  effects  of 
atmospheric  turbulence  are  investigated  [Ref,  6].  This 
model  uses  the  idea  of  a  long  term  optical  transfer  function 
(OTF)  when  considering  itmospneric  turbulence.  The  long 
term  OTF  results  from  simply  taking  an  image  of  sufficiently 
long  term  which  sees  effectively  all  possible  turbulence 
configurations. 

When  observing  taa  effect  of  diffraction  as  the  result 
of  the  finite  aperture  of  a  point  source,  it  can  be  seen 
that  the  resulting  image  is  not  uniform.  Considering  the 
point  spread  function  to  have  the  same  shipe  in  the  image 
plane  regardless  of  its  position,  the  image  function  is  the 
result  of   the  convolution   of  th=   source  function   and  the 
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optics  diffraction  function.    The   convolution  theorem,   as 
described  in  Fourier  transforn  theory,  yields 

i(v(x),v(y))=H  (v(x)  ,  v  (y)»  *o(v(x>  ,  v(y))  *M(v(x)  ,v  (y))    (2.1) 

where 

i(v(x),v(y))  =  Fourier  transforn  of  the  image  function 


H(v(x)  ,v  (y)  )  = 


Fourier 

function 


u  L  a  £1  a 


s for m  of  the  oo*:ic: 


o(v(x),v(y))  =  Fourier  transforn  of  the  object  function 


M(v(x)  f  (v(y))    = 


Modulation    Transfer    Function    of 
the    atmosphere 


v  =    spatial    frequency 

A  point  sourc=  such  as  a  Laser  can  be  analyzed  in  two 
dimensions  using  a  point  spreai  function.  Through  Fourier 
transform  theory,  the  image  point  spread  function  is 
transformed  into  a  two  dimensional  OTF  of  the  optical 
signal.  This  problem  may  be  simplified  by  scanning  the 
image  point  spread  function  using  a  vidicon  or  slit-scanning 
system.  In  applying  the  convolution  theorem,  the  Fourier 
transform  of  the  image  point  spread  function  is  multiplied 
with  the  Fourier  transform  of  the  optical  system  resulting 
in   the  Fourier      transform   of    the    overall      system.         The    Abel 
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transform  described  by  Griem#  is  applied  to  this  result  to 
re-transform  the  ois-diaensioaal  image  lsf  into  a  two 
dimensional    image   PSF    [Ref.    7]. 

As  demonstrated  by  Crittenden,  and  others,  a  numerical 
value  for  Cn2  may  be  obtained  by  curve  fitting  using  the 
following: 

M    =    exp(-21.49*Cn2*Z*f5/3*x-l/3)  (2.2) 

where 

H  =    MTF   of    the    atmosphere 

Z  =   Range   in  meters 

f  =   f*v   =   Angular  spatial   frequency   in   cycles/radian 

F   =    focal    length  of    the    op-ica.1    system 

v  =   linear    spatial   frequency    in    cycles/meter 

X  =    wavelength    in   meters 

Cn2    is   obtained   by   a   linear   regression   of 

ln(M)     =  -21.49*Cn2*z*fs/3*^-i/3  (2.3) 

where    Cn2    is  the   only    parameter. 
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III.   EXPERIMENTAL  PROCEDURE 

The  experiment  is  performs!  using  two  different  laser 
sources,  helium-r.eon  (He-Ne)  ani  gallium-arsenide  (Sa-As). 
Figure  3.1  shows  a  block  diagram  of  the  experimental  set-up 
used    for    both    lasers. 

The  measurement  of  Cn2  along  the  optical  path  is  made  by 
using  a  vidicon  ana  telescope  at  the  far  and  of  the  corridor 
in  the  basement  of  Spanagel  Hall.  Atmospheric  turbulence  is 
produced  by  nine  overhead  hot  air  ventilators.  The  optical 
equipment  similar  to  that  described  in  reference  3,  consists 
of  a  6  inch  diameter  CassegraLn  telescope  with  a  90  inch 
focal    length  [fief.    3]. 

Using  neutral  density  absorption  filters  to  attenuate 
the  intensity  of  the  laser  beam,  the  telescope  is 
illuminated  at  its  input  apertura.  The  light  image  from  the 
telescope  is  split  by  a  beam  splitter  with  one  beam  sent  to 
the  vidicon  and  the  other  to  a  slit  scanner  for  comparison 
of  the  two  systems.  The  light  is  then  transformed  to  an 
analog  signal  and  recorded  on  a  Panasonic (NV-9240)  video 
tape      recorder    for      processing.  The      data    from      the      slit 

scanning   system   are   recorded   on      a    precision    instrument    tape 
recorder   for      analysis.         After    recording   the      video    signal. 
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analysis  is  accomplished  using  a  Tektronix  468  Digital 
Storage  Oscilloscope  and  a  Hswlstt  Packard  9825  Computing 
System.  The  linearity  of  the  Panasonic  tape  recorder  was 
demonstrated  in  Crager's  thesis.  This  is  also  confirmed  by 
observing  the  images  as  real  time  and  recorded 
displays   [  Hef  -    3  ]. 

The  46  8  oscilloscope  is  set-up  in  the  following  manner. 
An  input  signal  to  channel  A  is  located  on  the  scope  using 
the    non-storage    mode.  The    signal    is    then      displayed    using 

the   A      INTENSITY   switch      of    the      horizontal    display.  This 

intensified  zone  is  used  to  position  the  3  sweep  (delayed) 
to  the  desired  location  within  the  A  sweep  interval  to 
obtain  an  expanded  view  of  a  waveform  for  examination.  Once 
the  waveform  is  centered  on  the  scope,  the  horizontal 
display  is  switched  to  B  DLY'D.  This  is  done  to  facilitate 
the  digital  storage  circuitry  time  base  by  using  the  setting 
of   the    B    TIME/DIV    switch. 

The  waveform  is  expanded  in  time  by  decreasing  the  A 
TIME/DIV  switch  setting  and  moving  the  waveform  back  to  the 
middle  of  the  scope  using  the  delay  time  position  control 
dial.  When  a  representative  waveform  is  obtained,  the  B 
TIME/DIV   switch    is      then  used    in    conjunction      with   the   delay 


time  position  control  dial  to  expand  the  time  zcal2  ,:",.~ii 
the  single  central  maximum  of  the  waveform  is  centered  and 
one  horizontal  sweep  is  displayed. 

Once  this  waveform  is  satisfactorily  obtained,  the  AVG 
storage  mode  is  selected.  The  468  will  average  the  input 
signal  for  a  selected  number  of  sweeps  and  display  the 
accumulated  waveform.  All  data  for  this  thesis  use  32 
sweeps  for  each  average.  The  463  is  now  ready  -c  transfer 
data  when  interrogate!  by  the  HP  9825. 

The  computer  is  the  controller  for  all  interfacing 
operations  with  only  a  minimal  amount  of  operator 
interaction.  The  operator  interface  is  mainly  to  ensure 
that  the  equipment  is  properly  set  up  and  to  select  if  data 
are  to  be  plotted.  Digitization  of  the  analog  signal  is 
accomplished  by  the  Tektronix  45"?.  Processed  waveform  data 
are  transferred  frnn  the  microprocessor  memory  to  the 
Storage  Display  RAM  [Sef.  8]. 

The  controlling  program  of  the  HP  9825  interrogates  the 
468  via  the  IEEE  488  interface  bus.  The  468  receives  the 
data  request  from  the  controller  and  sends  the  waveform 
message,  both  preamble  and  data.  The  waveform  message  is 
stored  in  the   calculator   memory   for  further   processing. 
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When  the  message  is  completed,  the  468  concludes  with  an  end 
of  instruction  tecmii ator  and  the  controller  takes  ccn-rol 
of  the  bus  again. 

Data  processing  o£  the  digitized  signal  begins  with  the 
data  being  stored  on  magnetic  tape  for  farther  use  by  the 
computer.  The  data  are  stored,  processed  and  plotted  by  th- 
main  program.  Subroutines  are  railed  as  necessary  for  their 
specific  uses.  The  two  signal  waveforms  are  recorded  and 
processed.  First,  the  signal  from  a  laser  beam  incident  at 
the  aperture  of  the  telescope  is  recorded.  Figures  3.2  and 
3.3  show  these  recorded  lata  waich  are  referred  to  as  the 
source  data.  Second,  for  calibration,  the  signal  from  a 
laser  beam  with  a  d if fraction  grating  in  place  at  the 
aperture  of  the  telescope  is  recorded.  The  entire  system  is 
calibrated  by  using  a  grating  that  consists  of  closely 
spaced  vertical  bars  in  front  of  the  telescope.  This 
produces  a  diffraction  pattern  in  the  image  plane  and  allows 
for  the  calculation  of  the  scale  factor.  These  recorded 
data  are  referred  to  as  the  scale  data  and  are  plotted  in 
Figures  3.4  and  3.5. 

During  the  recording  process,   the  signal  from  the  laser 
beam  is  recorded  on  video  tape   while  the  signal  is  observed 
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on  the  45  8  oscilloscope.  This  procedure  is  important  in 
data  taking  sines  lasar  alignment  and  telescope  focusing  are 
often  vary  sensitive  to  minor  movements  Df  either.  Upon 
completion  of  recording  the  lata  on  video  "ape,  the 
experiment  is  concerned  with  the  compilation  and  analysis  of 
data.  The  first  itea  to  be  determined  is  the  scale  factor. 
Once  the  scale  factor  has  been  da  terminal ,  it  will  remain 
constant  througnout  tae  calculations,  unless  the  telescope 
focal    length  is    changed. 

The  scale  factor  is  calculated  in  the  following  manner. 
The  number  of  points  between  tie  peaks  in  the  plot  of  the 
scale  data  (diffraction  grating  in  place  in  front  of  the 
telescope  aperture)  is  measured.  When  the  distance  between 
peaks,  the  spacing  between  the  bars  in  the  grating,  and  the 
wavelength  cf  the  laser  are  known,  the  scale  factor  can  be 
calculated    from    the   relation 

sin  (theta)    =  X /d  (3.  1) 

sin  (theta)    ~  theta  (3.2) 

(small   angle    approximation) 

sf    =    (theta) /ndp  (3.3) 

sf    =  X/  (d*nd?)  (3.4) 
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where 

sf   =   scale    factor   in   radians   per    point 

A       =    wavelength   in    asters 

d        =    spacing   between    lines   on   grating   in   meters 


ndp   =    number  of    points    between    oentral    maximum 
and   first    order    diffraction    peak 


This    standardizes   the    data    for   the      abscissa    of    the    plots   in 
radians    per   point. 

The  program  now  takes  the  transferred  data  and  computes 
a  point  spread  function  ae  shown  in  Figures  3.6  and  3.7. 
Figures  3.3  and  3.9  represent  the  point  spread  function 
after  integration  using  equation  (1.1)  to  obtain  a  line 
spread  function.  Next,  the  Fourier  transform  of  the  line 
spread  function  is  calculated.  These  curves  are  plotted  in 
Figures   3.10  and   3.11. 

The  diffraction  of  the  optios  is  computed  and  plotted  as 
shown  in  Figures  3.12  and  3.13.  Following  the  same  method 
as  before,  the  line  spread  function  of  the  optics  function 
is  calculated  and  pLotted  in  Figures  3.1'4  and  3.15.  The 
program  now  takes  the  Fourier  transform  of  the  optics 
function   LSF.         These      data    are    plotted   in      Figures    3.16    and 
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3.17.  The  Fourier  transforms  Df  the  system  and  optics  i;s 
divided  point  by  point  to  yieli  the  MTF  of  the  atmosphere, 
as      shown      in      Figures      3.18    anl      3.19.  The      MTF      of      the 

atmosphere    will    be   used  to      calculate   the   long-term    value   of 

Cn2  is  calculated  by  linear  regression  of  In (MTF)  versus 
f5/3.  The  slope  of  the  regression  is  proportional  to  Cn2. 
Two  simultaneous  equations  are  solved  using  Cramer's  rule 
and   the  matrix    ROM    of    the   HP    9825.       The    equations    used   are 

a*X    +    b*I    =    Y  (3.5) 

a*xz    ♦    b*X    =    Y*X  (3.6) 

where 

a   =    slope    of  the   curve 

b   =    intercept 

X  =    summation   of    (I*sf )  **(5/3) 

X**2   =   summation   of    (I  *sf)  **  (10/3) 

I   =    point    number    (total   number    =    256) 

Y  =   summation    of  natural   logarithm    of   MTF 

Y*X   =   summation   of    products   of    two    values 
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C  n2    is   obtained    from   the    icova    i  nr  ormat  ion    using    equation 


(2.3)  .      This  will   yield 


Cn2  =    a/(-21.  49*Z*f**(5/3)  *   **(-1/3)) 


(3.7) 


The  program  now  goes  to  the  prediction  phase  after 
calculating      C  n2      fsc      the        at  aosphere.  If      a      Gaussian 

distribution  for  an.  input  source  and  a  value  of  the  standard 
deviation  are  assunad  :r  known,  the  resulting  source 
function  can  be  calculate!.  A  plot  of  the  computed  source 
is  shown  in  Figure  3.20.  It  can  be  seen  from  this  plot,  that 
the  half- width  at  half-maximum  is  approximately  4 
microradian s.  As  before,  the  line  spread  function  of  the 
computed  source  is  calculated  aid  plotted.  This  is  shown  in 
Figure      3.2  1.  The      Fourier      transform      of      this      data      is 

calculated    and    plotted    in    Figure    3.22. 

The  program  now  multiplies  the  Fourier  transform  of  the 
computed  source  with  the  transform  of  the  system  (including 
the  atmosphere)  and  plots  the  result  in  Figures  3.23  and 
3.24.  Plots  of  the  inverse  Fourier  transform  are  shown  in 
Figures  3.2  5  and  3.25.  Next,  the  Abel  transform  is  computed 
and  plotted  in  Figures  3.27  and  3.23.  The  Abel  transform, 
as   described      in  Chapter      II,       transforms      a    one-dimensional 
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line  spread  function  into  a  tiro-dimensional  point  spread 
function.  Finally,  che  fraction  of  power  inside  a  circle  of 
radius  R  is  calculated  and  plotted  as  shown  in  Figures  3.29 
and   3-30.  Phis   is      the    fraotion    of      power    that      one   would 

expect  to  be  incident  on  a  target  using  the  measured  value 
of  atmospheric  turbulence  as  an  input  parameter.  A  detailed 
analysis  of  the  computer  progran  is  the  subject  of  the  next 
chapter. 
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IV.       COMPUTER    PROGRAM    IN    ALGORITHMIC    LANGUAGE 

The  program  preseited  hera  is  written  in  the  algorithmic 
language  as  outlined  by  Grahan  [Ref.  9].  It  is  intantied 
that  the  language  in  use  hara  be  machine  independent. 
Comments   appear   in    parentheses. 


Algorithm_Laser 

(Input    which    laser    is    usad.) 
Input   Laser^Waveiangth 
ESTRIUG<==    Type_of _Laser 

(Input  if  the  particular  lasar  source  data  are  recorded. 
The  following  lines  determine  if  the  laser  source  data  ara 
already  raccrded.  rf  the  souroe  data  ara  not  racorded  then 
the  program  goes  to  the  transfer  subroutine  to  allow 
interfacing  with  the  468  to  sffsct  a  data  transfar.  The 
parametars  passed  to  the  Plot  subroutine  have  the  following 
significance:  the  first  parameter  determines  the  horizontal 
index,  for  example  512  points  or  the  number  of  microradians. 
The  second  parameter  tells  tha  computer  which  file  tha  data 
are  recorded  on  so  taa  data  can  ba  loaded  into  memory.  The 
final    parameter     determines   the      labeling   of      the    horizontal 
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axis   either      in    data      points,      jiicrora  diar.s,         or    lir.i3      per 
microradian .) 

If   Laser_Source_no t_recoried    then 
Call    Transfer  (Source_File) 
Call    Plot  (512  ,Source_File,0) 
End   If 
(Input   if      the    particular      lassr    scale      data    are      recorded. 
Scale      data   are      recorded    with      the      diffraction   grating      in 
place-  If     the      scale      data    need      to      be      recorded,         the 

subroutine    Transfer    is    called   aid      the   scale    data   are   stored 
on    a    separate    file.) 

If   Lase r_Scale_not ^recorded    then 
Call   Transfer  (Scaie_Fiiei 
Call   Plot  (512  ,Scale_fiL  =  ,0) 
End    If 

Input  Scaie_Factor 
(Calculation  of  the  scale  factor  is  detailed  in  Chapter  III. 
Comparison  of  the  time  base  of  the  recorded  data  is 
necessary  to  ensure  that  the  lata  recorded  have  the  same 
time  scale  on  the  Dscilloscops.  If  the  time  scales  are 
different  the  program  stops  anl  displays  an  error  message. 
A   file   is   loaded,      the    waveform    preamble   is    searched    for    the 
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"XUUITS"  or  "XINCR"  osciilDScope  settings  and  the 
particular  data  are  stored.  This  is  done  for  each  of  tie  3 
recorded  files:  source,  scale,  and  collimated  data. 
Collimated  data  refer  to  the  same  laser  recorded  under 
controlled    conditions.  The    laser    beam      is    projected      in   a 

collimator  and  recorded  on  tape.  This  data,  could  be  used  to 
remove  any  non-linearities  of  the  video  recorder,  electronic 
equipment,  or  receiving  optics.  Data  elements  are  compared 
to   ensure   consistency    of   inputs.) 

Input   Laser_source 

3<==   po s(Laser_SDirce,  ' R:  ' ) 

R1<==    val( »  B+2'  ,3+U«) 

3<==    po s(Laser_SDurce,  'IT:  | 

ASTRING<==    Laser_5ource  (•  3  +  3  «  ,  •  3-m  ') 

Input  (Laser _Scaie» 

3<==   po s(Laser_Sca le, ' R: *) 

R2<==    val  (  •  B+2*  ,3+U') 

B<==    po s(Laser_Soale,  'IT:  ') 

3STRING<==    Laser_3cale  ( •  B+  3  •  ,  '  B+4  ') 

Input(Collimated_3  ource) 

B<==   pos(Collimated_So  urce, '  R:  •) 

R3<=  =    valCB+21  ,«  3+4) 

B<==    pos (Collin ate  d_SD  urce, ' IT: ') 
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CS?RING<==    Coll imated^So urea ( • B  +  3 ' ,  '3+4 ') 

If    R1  =  R2=R3      than 

If    ASTRING    =     BSTRING    =    ^STRING    then 

Input   RO  ATA(File_!Tumber) 
End  If 
Else 

Output     «ERR3R    IN    TIME    SCALE1 
Stop 
End    If 

A<==B<  =  =C<==0 
(Array   is    initialized    to   0) 
Do   for    I  =    1    to    51 2 

YD  ATA  (I)<  =  =    RDATA  (I) 
End   Do 
(This    section    takes    the    data    from    the    waveform    message    as 
transferred  from   the    468   and   removes   any    DC    background    from 
an    image   as  well  as    subtracting    the    minimum    value   in    order 
to    "zero"      the    data.     This    makes    y=0    as   the    starting    point 
for    both   calculation    and    plotting    purposes.    Data    up    to    51 
and  greater  than   463    are   summed    and   stored    in    B. 
This    represents   the      "wings"    of    the    curve.    This 
process   removes   the   DC   background    from   an    image.) 
Do   for    I  =    1    to    51 2 
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1<~    A  +  YDAIA(I) 
if   I   <   51    than 

B<==    YDfcTA(I)     +3 

End  If 

If  I   >   46  3   then 

3<==    YDATA(I)     +    B 
End  If 
End  'Do 
(Output   for   debugging) 

Ou-pur    A, 3 
(Average    of  the   sum    Df   data    values    in    the    "wings".) 

D<==   B/100 
(This   is   t he   total   sim    of    all    data    less    the    "wings"    value.) 

A<==    A-512*D 
(All    data    are    summed    and    store!    in    C.    Each    data 
is    compared   with   ths    previous    total    divided    by    2. 
This    will    find   the    ia  dex   value,    I,    for    which 
the   maximum  value    exists.) 
Do   for    I  =    1   t:    512 

YDATA(I)<  =  =    rDATA(I)-D 
C<=    YDATA  (I)     +    3 
If    C    <    A/2    then 
J<==    I 
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End  If 
End   Do 
(Output   for    debugging.) 
Output    C,J 
(The   curve   is    shifts!    so  that    ths    maximum    value    starts 
at   the   origin      of    ths    horizontal    axis.    This    is    called 
the   point    spread   function      and    is    rscordsi    as    RDATA.) 
Do    for    K   =     1    to    25  6 
IDATA  (K)<  =  = 
YDATA(  (J-K+  1)  mod512  +  1|  +YDATA  ((K*  J-1)  mod512*1) 
End   Do 
Do   for    I  =    257    to    512 

RDATA(I)<  =  =    3 
End    Do 
Do   for    I  =    1    to    256 

RDATA  (I)<  =  =    I  DMA  (I) 
End    Do 
(Determines  if    Tabular   or    piottsd    data    ars    desired.) 
If    PSF_OUIPUT_wa!ited   then 

Call    Table(RDATA> 
End    If 
If    PSF_PL3T_wantsd    thsn 

Call    Plot  (32,  RDArA,  1) 
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End   If 
(Calculates  the   line    spread   fuaction   from    the    point 
spread   function.) 

Call    LSF 
(Determines  if    tabuiir   or    plotted    data   are    desired.) 

If   LSF_OUTPUT_wanted   then 
Call    Table  (RD  ATA) 

Else   If   LSF_PLOT_wanted    thsa 
Call    Plot  (32,  RDAIA,  1) 

End   If 
(Calculation    of  Fourier  Transform    of    LSF.) 

Input    RDAIA 
(The    maximum    value    of    tha    array    IDATA    is    at    1=1.    Thi 
is   put   back  into   array   RDATA    and    folded  over    so   that 
IDATA  (257)    to    IDATA(5  12)       is    tie    mirror    image    of 
IDATA  (1)    to   IDATA(255).    Then    this    array   is    recorded 
as   RDATA   and   the  Fourier   transform   calculated.) 

Do   for    I  =    257    ta    512 
RDATA  (I)  <  =  =    3 

End   Do 

Do   for    I  =    1    to   256 

RDATA  (51U-I)  <==    RDATA(I) 

End   Do 
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HD  ATA  (2  57)  <==    R  DAI  A  (256) 
Call    FXFORM 

(Determines  if   tabular   or    plDttad    data   ars    desired.) 
If    FXFORM_OF_LSF_30TPJT_want3d    then 

Call    Table  (RD  ATA) 
Else   If   FXFOHM_OF_LSF_PLOr_wanted    then 

Call    Plot (32, RDArA,2) 
End  If 

(Dia  of    Cbscur/Dii    of   Dbj    Lens    in    meters, RATI0  =  . 064/.  16U) 
(This    is   the    measurement    of    the    ratio    of 
diameters  of   obscuration    to    telescope.) 
(R1=   Scale   of   data    in   microradians-l. 60) 
(Diameter  of   telescope   Obj    Lsns    in    meters ,  0BJECT= 
0. 1 64   meters) 
If    Laser_Source    =    He-Ne   then 
Wavelength<==    6.328*b-7 
Else 

Wavelength<=  =    9.35*e-7 
End   If 
(This    begins    the   calculation    o£    the    Fourisr    transform 
of   the   optics    using   the      Air/    Function.    All    of    the 
data    regarding   the    optics    are   recorded   on    a    separate    file.) 
(These   next  four   statements   fii    the    constants   in    the 
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argument   of  the   optics   function.) 
R1<=  =   5cale_Factor 
D<==    RATI0**2 
H<==    1-D 

Z<==   PI*Scale_Faot or*Dbject *  1 *e-6/Wave length 
Do   for    I  =    1    to    256 
Y<  =    Z*(I-1| 

If   Y   >    30   than 

RDATA  (II  <==    0 

Z  1  se 

RDATA (I|  <==     (AIRY(Y) -D*AIRY (Y*R ATIO) ) **2/H**2 

End  If 
K<=    I 
End  Do 

(Record    RDATA  on    separate   file) 
(Determines  if    a   plot    of    tha    optics   function    is    desirel.) 
If   OPTICS_FONCTIDN_PLOT    wanted    then 

Call    Plot  (32,  RDArA,  1) 
End  If 
(Calculates  the   line    spread    function    of   the    optics 
function. ) 

Call   LSF 
(Determines  if    a   plot    of    the    L3F    of    the   optics    function    is 
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desired.) 

If    ?LDT_OF_LSF_OPTICS   wantai    then 
Call    Plot  (32,  RDATA, 1) 
End   If 
Do   for    I  =    257    tD    512 

RDATA  (!)<==    3 
End    Do 
(Data   are    folded  over    again   anl    the    Fourier    transform 
calculated  .) 

Do   for    I  =    1    to    256 

RDATA  (514-1)  <==    RDATA(I) 
End    Do 

RDATA  (257)  <  =  =    RDATA  (256) 
(Record     RDATA) 

(Fourier   transform    of   the    LSF    of    Optics    is    calculated.) 
Call   FXFORM 

(Determines    if   plDt    of   Fourier    transform    of   Optics 
function    is    desirai .) 
If    PLOT_OF_FXFORM_OF_DPTICS    wanted  then 

Call    Plot (255 , RDATA, 2) 
End   If 

(Calculation    of    this    quotient    of    two    Fourier   transforms.) 
(This    section   takes    the    Fourier    transform    of    the   system 
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and   divides    it,    point    by    poiat,    by   the    Fourier    transform 

of    the   optics.    This    yields    tia    modulation    transfer 
function    (8TF)    of   the   itmospaere.) 

Input    RDATA_OF_FXFORM_0F_S0(JRCE 

Do   for    I  =    257    to    512 


RDATA  (!)<  =  =    D 


End    Do 

Do   for    I  =    1    to    256 

RDATA  (514-1)  C==    RDAIA(I) 
End   Do 

Input    RDArA_OF_FKFORM_OF_DPriCS 
Do   for    I  =    1    to    256 

RDATA(I)<  =  =    IDATA  (I) /RDATA  (I) 
End    Do 
(Determines  if    a   plot    of   the    HTF    is    desired.) 
If   PLDT_OF_MTF_OF_ATM3S_vantsd    then 

Call    Plot  (32,  RDATA,  1) 
End   If 

(Calculation   of   Cn2    by   curve    fitting. 
This    procedure   is    described    in    Chapter    III.) 
(All   "R"   variables    used    in   this    calculation   are 
initialized    to  0.) 
R17<==R  18<  =  =  R19<  =  =R20<==R2K=  =  R22<==R23<=  =  R2  4<==R25<==0 
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(Analysis    cf   data    shows   that   computar   tima    is   wasted 
b€yond   96   points.) 
Do   for    I  =    1    to    95 
(R18    -    when  summed    is    the     (f (ii  )  10/3.    This    is    the 
total   of   the    angular    spatial    frequencies    squared.) 
R1S<==     ((I-1|  *Scale_Fartor)  **(10/3) 
r  17<  =  =    R17    *■     R18 
(R19    -   when  summed    is    ths     (i(ii )  s/3.    This    is 
the  total    of   the  angular   spatiil    frequencies.) 
R19<==    ((I-1>  *Scale_Faotor) **  (5/3) 
R20<==    R19    +     R20 
(R21    -    is    the    total   number    of    points.) 
R2  1<==    I 

Do    While   RDM  A  (I)     >    0 
(R22   -    when  summed    is    the    total    of    the    product    of   the 
natural    logarithm    of    each    point    with    R19.) 
R22<==    in(RDATA(H  )  *R19 
R23<  =  =    S22    +    R23 
(R24    -   when  summed   is    the    total    of   the  natural    logarithm 
of   each   point. ) 

R24<  =  =    In  (RDATA(I)  ) 
R25<==    32U    ♦    R25 
End  Do 
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End    Do 
(Output   matrix    valuss    for    debugging.) 
Output     R17, R20, R21  ,R23, R25 
(Assign    value  to   matrix. 

This   section    takes    the    above    calculated    values    and 
sets   up    a  matrix.    The   matrix    equation    solved    is 
B    =    A"l*C.) 
<Cr~at=     3  aatzices    A  (2,  2)  ,3  (2  ,  1 )  ,Z  (2  , 1 ;  .  ) 


M1,1)<  = 
A  (1,2)  <== 
A  (2,1)  <  = 
A(2,2)  <  = 


R20 
R21 
R17 
R20 


INVERT, MAT RIX_A 
C(1,  1)  <=    R25 

C  (2,  1)  <  — R23 

B<==  A*c 
(Output  matrix  values  for  debugging.) 

Output  B(1  ,  1)  ,B  (1,  2) 
(Calculation  of  Cn2  using  equation  (2.  3)  .  ) 

R2  2<==  B(1,  1)/(-21  . 49* Range* R 20* (Wavelength**  (-. 3333  3) ) ) 

Output  'CNSQ=«,R22 

(The  program  now  starts  the  prediction  phase.   It 
calculates  a  source  assumed  to  have  a  Saussian 
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distribution    by    means   of   ths    computed   source    pattern 
below.    Using   similar  Fourier    transform    tschniques 
the   program    uses    the   calculated    value  of   Cn2    and    predicts 
the   power  incident    on   the    target.) 

(Computed  Source    Pa ttern : A=A3 *exp (-x**2/2*Sigma**2) 
(This   is   an    "arbitrary"    Gaussian    source   pattern 
with   a  standard   deviation  for    Sigma    set    equal    tc    2.) 
AC<==    &_ZERO 
C<==   Sigma*Scale_Factor 
Do    for    I  =    1   to    25  6 

F<  —     (1-1)  *5cale_Factor 
3<=    F**2/  (2*C**2) 
Do    While   3    >     13 

RDATA  (I)  <==    0 
End  Do 

RDATA  (I)<  =  =  A3*exp  (-G) 
End  Do 

(Record    RDATA.    This    records    the    computed    source.) 
(De-ermines    if   plot    of   computed    source   is    desired.) 
If   ?L0T_0F_C0MPUIED_S3URCE_wanted    then 

Call    Plot  (32,  RDATA,  0) 
End    If 
(Calculates  the   line    spread    function    of    the    computed 
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source- ) 

Call  LSF 
(Determines  if  a  plot  of  the  LSF  of  the  computed  source 
is  desired.) 

If  PL3T_OF_LSF_OF_COMP_SO(JR:E  =  1  then 

Call  Plot  (32,  RDAIA,0) 
End  If 

(Calculation  or  the  Fourier  c ransrorm  of  LSF  01  computed 
source. ) 
Do  for  I  =  257  -.o  512 

RD  ATA  (!)<  =  =  0 
2nd  Do 
Do  for  I  =  1  to  256 

RDATA  (51U-I)  <==  RDAIA(I) 
End  Do 

RDATA(257<==    RDAIA(256) 
Call    FXFORM 
(Determines  if   a   plot    of   the    Fourier   transform    of   the 
LSF    of   the   computed   source    is    desired.) 

If    PLOT_0F_FXF0Rtf_OFj:OMParSD_SOURCE_warited    then 

Call    Plot  (255  , RDATA, 2) 
End    If 
(Calculation    of   the    product    of    two    Fourier    transforms.) 
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(This   is   the    Fourier   transform    of   the   computed    scu:c3 
with   the    MTF   cf    ttis    atmosphsc a. ) 
Do    for    I  =    1    to    25  6 
RDATA(I)<=  = 

RDATA  (I)_C0SPUTED_S3aRCE*RDATA  (I)  _FX  FK_0FJ1TF 
End    Do 

(Product   of    FXF0RM5    of    Sourca* Atmosphera*0 ?t ics. ) 
RDATA (!)<== 

RDATA(  1)_?RODUCT_D  F_2_FXFMS* RDATA  (1)  _FXFM_OF_OPTICS 
Do    for    T  =      2   to    2  56 
RDATA  (I)<=  = 

2D  ATA  (I)  _?RD0UCT_0F_2_?XFKS*RDATA  (I)  _FXF**_OF_0?TICS 
RDATA  (514-1)  <==RDATA(I) 
End  Do 

RDATA(257)  <==    RDATA(256) 
(The    above   data   predicts    what    ths    Fourier    transform 
of    the   entire    system    isf    using    the    calculated 
value   of   Cn2    for   the    atmospheric    turbulence.) 
(Determines  if   a   plot    of   ths    rasult   of    Fourier    transform 
is   desirad.) 

If    P  LOT  J3F_FXF0RM_  PROD  OCT  Sjf  anted   then 

Call    Plot  (32,  RDATA, 2| 
End   If 
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(Inverse   Fourier    transform    gives    target    LSF.) 
Call    INVERSE_FXFDRM 

(Determines  if    a   plat    of   the    inverse    Fourier    transform 
is    desired.) 

If    ?L0T_OF_INVERSE_FXFROM_wanted   then 

Call    Plot(32,  RDArA,0> 
End  If 
(Converts    target's    Die   dimensional   LSF   to    a 
two    dimensional  PSF    by   Abel   transform.) 

Call    ABEL 
(Determines  if   a   plot    of   the   Abel    transfoca    is    desired.) 
If    PL0T_0F_A5EL_XF0RM_wanted    then 

:all    Plot  (32,  RDAIA,0| 
End  If 

(Calculates    the    fraction   of   power    inside    circle   of   radius 
R.    This    predicts    the    fraction    of    power    that    will    be 
incident    on    target.  ) 
RDATA(1)<  =  =    0.25*PI*RDATA(1> 
Do    for    I  =    2   to    256 

R DATA  (!)<  =  =    2*PI*RDAIMI)  +RDATA  (1-1) 
End   Do 
(Determines  if   a   plot    of   the    fraction   of   power    is   desired.) 
If    PLOT_OF_POWER_* anted    thea 
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Call    ?lot(32,  RDArA,Ql 
End   If 
End    Laser 

(The      following      subroutines      i"£    used        in      the      Algorithm 
Laser.) 

Subroutine_Transf Br (File_ Number) 

(The   transfer    subroutine   gets    raw    data,    both    Draamble 
and    data    from    the    waveform    massage    sent    by    the 
U6S.    It    processes    the    message    by    finding    the    minimum 
value    of   the    array    then   subtract ing • this    value    from 
each   element    in  the   array.    This    "zsros"    tha    array.) 
Output    'Ensure   E:j:iipm=nt    set    up    properly' 
Output    'Continue    *  hen    ready' 
E<==    (-1) 
Do    Whila  S    =   -1 
("DATA"    is    the   I/O    buffer   wher=    the    data    from   the   468    are 
sent.    The    status   of    the   buffer    is    read   while    data    are    being 
transferred,    tfhen   tha    transfar    is    complete,       the    interface 
is   cleared.) 

Data<==   Trans ferred_Data 

If   Data_Traas fer_Compl2t e    then 

E<==    0 
End  If 
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Clear_ Inter  face 
N<==    1 

(This    is    where   heaier  format    stops    and   data    begin.) 
A<==    pos(Data,  •  %»  > 
Do   for    I  =    A+2    to    687    by    16 
Do   for   J   =    1    to    16 

If   J+I>5  88    then 
RJ<==    3 
N<  =  =    N+1  ' 
End   If 


RJ<=  =    nam(Data  (I*  J)  ) 
(The   numerical    valu25    of    each    aiement    are    stored    in    YD\TA.) 
YDAIA  (M|   <==    RJ 
N<==    Nf1 
End    Do 
End   Do 

E<==    minimum  (YD Ar A) 
Do   for    I  =    1    to   51 2 

YDATA  (I)  <==    YDMA(I)-E- 
RDATA  (I)  <==    YDAT  A(I| 
End   Do 
Store_RDATA (File_N umber) 


73 


Return 


End   Transfer 


(This   subroutine  calculates   th  =    Line   Spread    Function 
for   the   previously   recorded   point    spread   function 
by    using    equation    (1.1).) 

Subrout  ine_LSF 
(Depending    on    which    file    number    has    been    passed    down 
from    the   calling   subroutine,    either    scale    or    source 
data   are    loaded.) 

If    File    Number    =   Scale    Data    then 


Input    Scale   Data 


RD  ATA<==SCALE     DAr  A 


Else 

Input    RDATA 
End   If 
Do   for    I  =    1    to    51 2 

ID  ATA  (I)  <==    3 
End    Do 
Do   for    I  =    1    to    25  6 

IDATA(I)<  =  =    RDATA  (I) 
End   Do 
(Plots  have  shown    that   computer    time    is    Ids-    and    no 
valuable   information    is   gained      beyond   about    24    points.) 


7  it 


Do    for    I   =    1    to    21 


J<==    1 


Output   I 
Q<=    IDATA  (I) 
Do   While   a    <    24 

R<==    S2RT  (I*I  +  J*J» 

Q  +  2*(  (1-  fraction  (R)  )  *I  (int  (?)  )  ^fraction  (?)  *I  (int  (E)  +  1)  ) 
J<==   J+1 

End  Do 

IDATA(I)<  =  =    2 
End   Do 
Do    for    I  =    1    to    2'4 

RD  ATA  (I)  <  =  =    IDATi  (I) 
End    Do 
If   File_Number   =   Scale^Da^a    then 

■     St ore_Scaia_D ata 
Else 

Store_RDATA 
End   If 
Return 
End   LSF 
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(This   subroutine   calculates   ths    Fourier    transform 
of    the   given    data    using   the   CDoley-Tukey    Algorithm, 
If   the   inverse   statement    is   true,    then   ths    inverse 
Fourier    transform    is    calculate!.) 
Subroutine_FXFORM 

Set_Rad ian^Mode 
(2**9=512  which    is    the   number    of   points.) 

N<==   9 

Do   for    I  =    1    to    51 2 
IDATA  (I)<  =  =    3 

End   Do 

If    File    Number    =    Scale    Data    then 


Input    Scale_3ata 


Else 


Input    RDAIA 
End   If 
T<==    PI/2**  (N-1) 

Do   for  J   =    3    to    2**  (3-1)  -1 
(BI   is   a   bit   inversion   subroutine.) 
Call   31 (J,P, N-1) 
c<==  cds  (P*r) 

If   INVERSE   FXFORM    =    True   then 


Flg7<==    1 
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Fig7<==    0 
End   If 

P<==    sia  (P*D  *(1-2*Flg7) 
Do    for    I    =    2*R0*J*-1    to    2*R0*J+R0 
R1<==    RDArA(I) 
R2<==    RDATA(H-RO) 
R3<==    IDATA  (I) 
R'4<«    IDATA  (I  +  RO) 
RDATA(I)<  =  =    31*R2*OB4*P 
IDiVTA(I)<  =  =    R3+R4*C-R2*P 
RD&TA  (I  +  RO)  <=  =    B1-B2*C-Rt**P 
IDATA(I  +  R0|  <  =  =    R3-R4*OR2*P 
End   Do 
End  Do 
Output    H 
End   Do 

(This   s  act  ion  is   for   re-ordsring   the   block.) 
Do   for    I  =    0   to    2**N-1 
Call    BI(I,J,tf) 
Do   While   I-J    <   or    =    D 
Do   Whils    I  #    J 


?<==    RDATA(I>1)  /SQRT(2**N) 
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Z<==    IDATA  (I*  1)  /SQRT  (2**N) 
RDATA(I-H)  <  =  =    RDATA(J  +  1| 
ID1TA(J+1)<"    IDATA(J+1) 
RDATA  (J  +  1)  <  =  =    P 
IDATA  (J+1)  <  =  =    Z 
End   Do 

RDATA  (If  1)  <==    RDATA  (1+1)  /SQRT  (2**N) 
IDATA  (l>  1)  <=  =    IDATA  (H-1)/SQRI  (2**N) 
E nd    Do 

If  File_N umber   =    Scal?_Da-a    then 
Stors_S::ale_Data  (Fils_Number) 
Else 

Stor=_RD  &TA( Filers umber) 
End   If 
End   Dd 

Set_Degrees_Mode 
Return 
End    FXFORM 

(This   subroutine   takas   the    binary    number   ?1    containing 
P3    bits   and  inverts    it    end   for    ?nd,    e.g.    310111    becomes 
111010.) 

Subroutine   BI(P1,?2,P3) 
P2<==    0 
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PU<  =  =    p  1 

Do   for    2  =    1    t:    ?3 

P4<==    ?4/2 

?2<==    2*P2 

If   fraction(P4)    »    0   than 
P2<==    ?2    ♦    1 

End   If 

P4<==    int(P'4| 
End    Do 
Return 
End    31 


(This   function    calculates   Airy  ( x)  =2*J1  (x)  /x    where   J1  (x) 
is   the    Basse!    function   of    ordac    one.) 
Function   AIRY(P1> 
If   P1    <   0    then 

Output     'ERRDR -ARGUMENT    LESS    THAN    3« 
Step 
End  If 

If   PI    ■  0    then 
R  4  <==     1 
Return    RU 
End   If 
R5<==   0 
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If   P1    >    15    then 

R6<==    90    *    P1/2 

5 1 2<==    1.4*P1     *    60/P1 

End   If 

If    P1    <   5    then 

R6<=  =    20    *■    10*Pl-P1**2/3 


R12<  =  =    6    ♦    P1 


Else 


R6<==    20    +    13*P1-P1**2/3 
R1  2<==    1.4*P1     ♦    50/P1 
End   If 

R12<==    maximum  (lit  (P.  12)  ,int  (3+P1/4) ) 
Do    for    a   =    R12    to    R6    by    3 
R8<==    1*e-23 
R13<==R14<  =  =    0 
If   M/2    =   int  (M/2)     than 

Flg10<==    0 
Else 

Flg10<==    1 
End  If 
Do   for   J  =    1    to   S-2 

R15<==    2*(M-J) *R9/P1-R13 
R13<==    R8 
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R8<==S15 

If   M-J-2    =   0    than 

R'4<==    B15 
End   If 
If   Flg10    =   0    than 

Flg10<==    1 

ElS9 

Fig  10<==    0 

End   If 


R14<==    R  14  +  2*R8*FL?10 


End  Do 

R1  5<  =  =    2*R8/?  1-R1  3 

R1 4<==    R14+R1 5 

R4<==R4/R1U 

If    (abs  (RU-R5)  -abs(R'4*1  *?-6)  )     <    or    =    0    then 
RU<==    2*RU/P1 

End  If 

Return    R1 

R5<==    R4 
End    Do 

Output     'ACCURACY    NOT    OBTAINED1 
Return    RU 


End    AIRY 
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(The   subroutine   Abel    takss   a    019-dimensional 
line    spread  function    and   calculates   a   two 
dimensional  point   spread   function.) 
Subroutine    ABEL  (Fila_ Number) 
N<==    RDATA  (1) 

RDATA( 1)<==    1.4*RD ATA( 1) -1 . 3*RDATA (2) ♦. 4* RDATA (3) 
Do    for    I  =    2    to    5 '4 
M<=    RDATA  (I> 

RD  ATA  (I)  <  =  =    .  4*N+.2*H-.  6*BDATA(H-1) 
N<=    M 
End   Do 
Do   for    I  =    1   to    5 '4 

RD  ATA  (I)  <  =  =    RDATA  (I)  /(2*SQRT  (  (I  +  .  1)  **2-I*I)  ) 
Do    for    J    =    1     to    5  4 
RDATA  (I)  <=  = 
RDATA  (I)  +RDATA  (J) /SQRT  (  ( J* .  1 )  *  *2  - 1  *I) 
End   Do 

RDATA  (I)<==    RDATA  (I) /PI 
Output    I 
End   Do 
Return 
End   ABEL 
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(This   subroutine  prints   out   data    in    tabular    form.) 
Subroutine_TABLE(Fila_Nuaber) 

Input   RDATA  (File_tf  umber) 

Do    for    I  =    1    to    32 

Do   for   J   =    1    to    15 
Output    RDATM  16  (1-1) +J» 

End    Do 

Output     RDATA(16*I| 

Return 
End   TABLE 

(This    subroutine   plDts   the   desired    data.) 

(Each   time   a    plot    is    called,    tils    subroutine 

plots   the    particular    graph   and    labels   it. 

The   horizontal    axis    is    labelad    according   to    the    first 

and  last    parameters   passed   by    the    calling    subroutine.) 

Subroutine_PLOT (P1 ,Fi le_N umber , ?3) 

Input    RDATA (File_tJ  umber) 

A<==    0 

3<==    P1 

C<==    minimum  (RDAI&) 

D<==    maximum  (RDArA) 

XMIN<==    -.  1*  (B-A) 
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XMAX<==    B«-.05*(3-fc) 
YMIN<==    C-.  1*  (D-Z) 
YMAX<==    D+. 05*(D-3) 
E<==    B 
F<==    10 
(Using      the    valua      of   P1  ,        ttia      horizontal    increment      foi 
plotting   and  labeling    is   determined.) 
If   ?1    =    512    than 

3<==    6  4 
Else   If    P1    =    256    then 

3<==    3  2 
Else   If    ?1    =64    then 

3<  =    8 
Else 

3<==    4 
End   If 
Plot    B,C,1 
(This    lines   and  places   tic    macks   on    the    horizontal    axis.) 
Do   for    I  =    E  to    0    by   -G 

Plot    I*  (B-A>  /E+A,  C,  2 

Plot    I*  (B-A)  /E+A,C+  (D-3)  /150,2 

Plot    I*(B-A|  /E+A,C,2 
End    Do 
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(This   lines   and  piires  tic   aar'vs   en    the   vertical   axis.) 
Do   for    I  =    0  to   ? 

Plot    A,I*(D-C)/F+C,2 
Plot    A  +  (B-A)/150,I  *(D-C)  /F  +  Z,  2 
Plot    k,  I*(D-C)  /F«-:  ,2 
End    Do 


Charact  er_Size    1.2,1, .7,0 


rtical  axis. ) 


(This  numerically  libels  the  n 
Do  for  I  =  F  to  3  by  -1 

Plot  A-.075*  (B-A)  ,1*  (D-C)  /F  +  C,1 

Label  I/F 

End  Dc 
(This  sets  up  a  character  string  for  labeling  the  vertical 


axis. ) 


If  P3  =  0  then 


ASTRING<=  =    'DATA    POINTS" 


L<=    1 


Else   If    P3    =   2    then 

aSTRING<~    ' LINES/MIZRDRADIAN* 

L<~    1/(2*L| 
End   If 
(This   numerically    libels   tha    horizontal   axis.) 
Do    for    I  =    0    to    E    by    3 
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Plot    A*  (I/E-.  025)  *  (B-Al  ,  C-  .  0  25*  (D-3)  ,  1 

Label    I*L 
End   Do 
If   P3    #    0    then 

L<  —    Scale_Fictor 
End  If 
(This    labels    the    horizontal    axis.) 
Plot    .  4*(B-A)+A,.  0  5*  (D-C)  *C,  1 
Label    A  STRING 

Plot    -.  07*  (B-A)  +A,  .3*(D-C)  +  C,  1 
Characxer_Size    =    1.2,1,-7,90 
A3TRING<==     'NORMALIZED    INTENSITY" 
Label    ASTRING 

Charact er_Size    =   .5,1,1.5,0 
(This    plots    data.) 
I<==    0 

Plot   I,  maximum (RD ATA  (I)  )  ,1 
Do   for    I  =    1   to   P1 

Plot    I-1,RDArA(I) 
End   Do 
If   P3    =   0    then 

Input    PLOT  _L  A  BEL 

ISTRING<==    PLOT_LABEL 
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End    If 

Plot    .6*B,.9*D,1 
Label    • PLOT    OF' 
Plot    .  6*B, . 87*Df 1 
Label    ISTRING 
Plot    .  6  *B, .  8<4*D, 1 
Label    ESTRINGr    •    Laser' 
Return 
End    PLOT 
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V.   CONCLUSION 

The  work  reports!  in  this  thesis  supports  the  uodel 
predicted  and  measured  in  Crittenden,  and  others,  for  the 
long  exposure  case  [Ref.  2  ]•  The  vidicon,  in  replacing  the 
mechanical  slit  scanning  system,  shows  hd  degrading  of  the 
signal  data.  It  also  demonstrates  that  a  gDod  approximation 
for  the  point  spread  function  nay  be  made  by  recording  a 
single  TV  line  through  a  laser  spot.  This  line  is  then  used 
to  calculate  the  one-i imensional  line  spread  function.  The 
linearity  of  the  video  tape  recorder  is  seen  in  the  results 
of  the  MTF  and  Cn2.  Tha  463  Dsoillcscops  is  the  workhorse 
for  tha  entire  systen.  It  affectively  displays,  stores, 
digitizes,  and  transfers  data  on  a  real  time  basis. 

The  overall  systen  does  not  have  the  capability  of  the 
system  described  by  Cragec.  However,  the  relative 
simplicity  cf  the  structured  programming  technique  coupled 
with  a  lower  eguipnent  cost  demonstrates  that  comparable 
measurements  and  data  evaluation  using  this  equipment  can  be 
made.  Further,  investigation  should  include  use  of  the  Data 
Precision  6  000  digital  waveform  analyzer  thereby  allowing 
data  to  be  sampled  and  analyzed  on  a  near  real  time 
basis. 
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APPENDIX  A 


0:     "HE-NE:  SOURCE   ON    9  ,  SCALE   ON    7,&    COLL  I  MATED    ON    5": 
1:    "GA-AS.-SOURCE   ON    8  , SCALE    ON    6 .&CCLLIMATED    ON    4": 
2:    ent    "SELECT   CODE    FOR    PRINTER" , A;dev    "print"  ,A 
3:    fmt    l,z  , c  ;tmt    2,  td  .1  ,z;  fmt    3,18.1 
4:    dim   I$.[1024]  ,R[512]  ,Y  [  512]  ,C$  [20]  ,D$  [20]  ,E$[2Q] 
5:    aim    I  [512]  ,A$.[32]  ,B$.[20]  ,A[2,2]  ,fi(2,l]  ,C[2,1] 
6:    but    "DATA"  ,I$,.3;ina    X  ,1  ,R;  0*R*W+V-KJ;  1+L;  but    "I"rI$r3 
7:    beep; asp    "EN1EF      LA SEE    USED"; wait    150  0 
8:    ent    "HE-NE=1 ,GA-AS=0 " ,w 
9:     it    W=1;"HE-NE"*>E$;9+C;  jmp    2 
10:     "GA-AS"-E$;  8-Q 

11:    fceep;dsp    £$&"    SOURCE    DATA    RECORDED?"  ; wai t    1500 
12:    ent    "1  =  YES  ,  0=NO"  ,V 

13:     it    V=0   ana    W=l;  stg    0 ; ell  "'TRANSFER  '  (Q)  ;  jmp    2 
14:     if    V.=0    and    to=0;£tg    0  ;cll  "  'TRANSFER  '  (Q) 
15:    beep;ent    "PLOT    OF    SOURCE? , 1=YES  ,G=NC"  ,R 
16:     if    R=l;cli  "PLOT" (512,C,G) 

17:    beep;asp    E$&"    SCALE    DATA    RECORDED? "; wa it    1500 
18:    beep;ent    "1=YES , 0=NO" , U 

18:    ir    U=0   ana    w=l;stg    0 ;  ell  "'TRANSFER '  (7+OJ  ;  jmp    2 
20:    if    U=G    ana    Vv  =  0;  sfg    0; ell  "'TRANSFER ' (6+G) 
21:     aeep;asp    "WHAT    IS    SCALE    FACTOR?"  ;  Vv ait    1000 
22:    ent    "SCALE    FACTOR=? ; 0    GETS    PLOT" ,S 
23:    it    S=G    ana    Q#0; Q-2*Q;cll  '  'PLOT  '  (512 ,Q, Q) 
24:     if    S  =  0    ana    U  =  0 ; ell  '  'PLOT  '( 512 ,Q , 0) 
25:    it    S=0;3mp  .-3 
26:    S-L 

27:    if    W=l;  9+.T+Q;* jmp    2 
28  :    8+.T+C. 
29:     laf    T,I$,R[*] 
.30:    pos<I$, "R:") -B 
.31:    vai(I$[B+2,E+4] )*rl 
.32:    pos(I$,  "IT:".) +6 
33:    IS[B+.3,B+4]-A$ 
34:    lot  .T-2,I$,R[*] 
35:    pcs(I$,"R:")+B 
.26:    val(I$[B  +  2,B+4].)+r2 
37:    pos(I$,"LT:'\)+B 
38:    I$[B+3,B+4]+B$ 
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39:    lcir    T-2  ,I$,R[*J 
40:    pos(I$,"R:" )  -B 
41:    val(I$[B+2,B+4] )+r3 
42:    pos(I$/"LT:,\)+B 

43:    I$[B+3,B+4]-C$ 

44:    if    r2  =  rl    ana   A&=B$    ana    r2=r3    ana   B$.=C$;jirp    2 

45:    beep; asp    "ERROR    IN   TIKE   SCALE " ; stp 

46:     "B":ldf    Q,  1$  ,R[*  ] ;  ina    Y;ara   R-Y;C-A-B-k: 

47:    for    1=1    tc    512 

4b:    A+Y[I]-A;if    K51;B+Y[I]*B 

49:     if    I>463;B+Y["l]+B 

50  :    next    I ;  pr  t   A  ,B 

51:    B/100-»»D;A-512*D-A 

52:    tor    1  =  1    tc   512;  Y  [  I  ]  -D  +  Y  [  I J  ;  C  +  Y :[  I J  *C 

53:     if   C<A/2;  I  *  J 

54:    next    I;prt    C,J 

55:     for    K=l    to    256 

56:     Y[  (J.-K+L)ir.ca512  +  1]  +Y[  (K+J-l)  ir.cc512  +  l]  *I[K] 

57:    next    K;ma    R;ara    I-*R 

58:     "IMAGE    POIN.T    SPREAD    FCN"-I$;rcr     10,I$,R[*] 

59:    beep;ent    "PSF    CJTPUT=1   AND/OR    CCM"  ,Z 

60:     beep;ent    "PSF    PLOT  =  l   AND/OR    CONT"  ,Y 

61:    if    Z=l;  ell  ['.TABLE  '(10) 

62:     if    Y  =  l;cli  "  'PLOT ' ( 32 ,10 , 1 ) 

63:    ir,in(l  [*]  )*r22;max  (I[*].)+r23 

64:    prt    r22,r23 

65:    oeep;ent    "LSF    OUTPUT.  =  l  AND/OR    CONT"  ,Z 

66:  beep; ent  "LSF  PLOT=l  AND/OR  CONT"  ,Y 

67:  ell  "LSF '; "IMAGE  LINE  SPREAD  FCN"-I$;rcf  10,I$,R[*] 

68:  it  Z  =  l;cll  ".TABLE  '(10) 

69:  if  Y  =  l;cll  "PLCT"(32  ,10,1) 

70:  "CALCULATION  OF  FXFORM  OF  LSF":laf  10 , 1  $ ,R [ *]  ; ina  I 

71:  ara  R+I 

72:    I.[l]-R[l]  ;for    1  =  2    tc    256 

7  3:    I [I] +R[I]*R[ 514-1] ;next    I 

74:    R[256]*R[257] ;rct    10,I$,R[*] 

75:    Oil  "FXFORM";"  FXFORM   OF    IMAGE    LSF"+I$;rcf    10,I$,R[* 

76:     beep; ent    "FXFORM   OF   LSF    OQTPUT.  =  l   AND/CR    CONT"  ,Z 

77:    beep;ent    "FXFORM   OF    LSF    PLOT  =  l   AND/OR    CCNT"  ,Y 

78:    if    Z  =  l;cll  "TABLE"  (10) 

79:    if   Y=l;cll  "PLOT" (32 , 10 , 2) 

80:     "B=DIA    OF    CBSCUR/DIA   OF    OBJ    LENS    IN    METERS": 

81:     "B=0.G64/.164": 

82:     "rl=SCALE    OF    DATA    IN   MICRORADIANS , rl+1 . 6  0" : 

83:     "0=DIA   OF    OBJ    LENS    IN   METERS ;O=0. 1 64M" : 

84:     "K=toAVELENGTH    IN    METERS": 
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85:     "CALC    CF    DIFFRACTION    LIMIT    POINT    SPREAD    FCN" : 

bo:    ina    R;i£    w«l?  6.328e-7+fo;  jmp    2 

87:    9.05e-7-»-W 

68:     . 064/. 164-E; .164*0 ;L-rl 

89:    B  *2+D  ;  1-D+H;  7i*rl*0*le- 6/to-Z 

90:    for    L=l    tc    256 

91:    Z*(I-l)+Y;i£    Y>50;^tc    +3 

92:     ("AIRY"  (Y)-D*"AIRY "  (Y*B)  )  ~2/H~2-R[I] 

93:    I+K;gto   +2 

94:    0+R[I] 

95:     txa  .5;Qsp    R[I];next    I;"OPTICS    FUNCTION "+I $ 

96:    ret    11,I$,R[*] 

97:  beep;ent  "OPTICS  FCN  PLOT=l  AND/OR  CONT"  ,Y 

98:  it  X«l;cll  "'PLOT* '(32,11,1) 

99:  stg  5 ;cll  "LSF"; "LSF  OF  OPTICS  FCN " -I  $ ; ret  11,I$,R[*] 

100:  beep;ent  "PLOT  OF  LSF  OF  OPTICS  FCN=1  AND/CR  CGNT",Y 

101:    it    Y=l;cll  "PLOT  '  (32 ,11, 1) 

102:    stg    5;  ina    I;  ara    R+I;  I.[  1]  *R[  1] 

103:    tor    1=2    tc    256 ; I [ I ] *R [I ] -R[ 514-1] ; next    I 

104:    R[256]  +R[257]  ;rcr    I1,I$,R[*] 

105:     ell  "FXFCRM"; "FXFCRM    CF    LSF    OF    OPTICS" +1$; rxf    11,I$,R[*; 

106:    ceep;ent    "PLOT   CF    FXFM  OF   OPTICS=l  AND/OR    CONI",Y 

107:    ir    Y«l;cll  "PLOT  '(256,11,2) 

106:     "CALCULATION    OF    QUOTIENT   OF    TtoG    FOURIER    TRANSFORMS": 

109:     laf    10 ,I$,R[*J ; ina    I; ara   R-I 

110:     idf    11,I$,R[*] 

111:    tor    1=1    tc    25o 

112:     I [I]/R[L] -R[I]  ;next    I 

113:     "MTF    CF    SYSTEM" +1$ ; ret    10,I$,R[*] 

114:    beepjent    "PLOT    OF    MTF    OF    SYS=1   AND/CR    CONT" ,Y 

115:    it    Y=l;cll  "PLOT  '  (32, 10 , 1) 

116:    0-*rl7-»rI8-*rl9->r20-*r21-r2  2-r2  3-r2  4-r2  5 

117:    tor    1=1    tc   96 

118:     ( (I-l)*L)~(10/3)-rl8 

119:    rl8+rl7+rl7 

120:     ( (I-1)*L) "(5/3) +rl9 

121:    rl9+r20*r20 

122:    I-r21 

123:    if    R[I]<=0;gtc    "CC" 

124:     (In (R[I] ) *r24) *rl9+r23-r23 

125:    r24+r25-r25 

126:     "CC":next    I 

127:    prt   rl7,r  20  ,r2*l,r23  ,r25 

128:    r 20 -A [1,1] ;r 21*A [ 1  ,2] ; rl7-A[2  ,1]  ; r2 0+A[ 2, 2 ] 

129:    r25*C[l,l]  ;r23-Kl[2,l] 

130:    inv   A+A;ir.at   A*C+B;£lt    5;aprt   B 

131:     "CALCULATION    OF    CNSQ": 

132:    E[l,l]/  (-21. 49*145*  r20*W~  (.-.33333)  )  *r2 2 

133:    prt    "CNSQ-"  ,r22; stp 
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134:  cfg  ; laf  2 

135:  "TRANSFER" : beep; asp  "ENSURE  EQUIP.  SET  UP  PROPERLY" 

136:  wait  150  0 

137:  asp  "PRESS  CONTINUE  WHEN  READY" ; Stp 

138:  dsp  "waiting" ; wait  1500 

139:  our  "DATA"; 0*2 

140:  tfr  703 /'DATA" 

141:  asp  "transfer ring" 

142:  rds("DATA")  -E ;  it  E=-l;]ir.p  0 

143:  clr  703 

144:  cisp    "setting    output" 

145:  pos (1$, "%")*A;1*N 

146:  for    I=A  +  2    tc    687    by    16 

147:  ror   J=l    tc   16;if   J+I>688; 0*rJ ; jmp   3 

148:  num|l$[I+J])*rJ 

149:  rJ+Y[N] 

150:  N+1*N 

151:  next   J 

152:  next    I 

153:  min(Y[*] ) -E;ror    1=1    tc    512; Y[I ] -E+Y [I ] ; next    I 

154:  ina    R;ara    Y  +  h;  rci    pi , I  $  ,R [ *]  ; r et 

155:  "LSF":if    tlg5;iaf    11 ,1$  ,R[*j  ;  jmp    2 

156:  laf    10,1$, R[*] 

157:  ina    I;ara   R+I;for    1=1    tc    24 

158:  1-J;csp    I ; L [ I ] -C 

159:  vd*I+0*J)*R 

16  0:  2*  ( (1-frc  (R)  )*I[mt(R)  J+frc(R)*I[int(R)+lj)4g*C 

161:  J+1+J;if    R<24;jirlp   -2 

162:  Q*I[I] 

163:  next    I ; ina    R;ara    I*R; if    rlg5;rcf    11, 1$  ,R[*J  ;c£g     ;ret 

164:  ret   pi 

165:  "FXFORM" :rad;9*N; ina    I ; if    rlg5;laf    11,1$ ,R[*] ; jmp    2 

166:  laf    10,I$,R[*] 

16  7:  TT/2*  (N-l)  *T 

168:  for    M=l    tc    N  ;  2~  (K-k)  -rO 

169:  for    J  =  0    tc    2"  (M-l).-l;  ell  '  'BI  '( J,P  tN-l) 

17  0:  ccs  (P*T)  *C  ;sm(P*T)  *  (1-2* tig  7)  *P 
171:  tor    I=2*rO*J+l    tc   2*rG*J+rG 
172:  R[I]+rl;R[I+rU]*r2 

173:  I[I]*r3;I[I+rO]*r4 

174:  rl+r 2*C+r4*P-R[ I] ; r 5+r4*C-r2*P*I [ I] 

17  5:  rl-r2*C-r4*P-R[I+rC] ;  r3-r  4*C  +  r2*P-*I  [  I  +  r0  ] 

176:  next  I; next  J; asp  M; next  M 

177:  for  1=0  to  2*N-l;cll  "'BI  '( I , J ,N) 

178:  ir  I-J>0;gto  "SB" 

179:  if  I=J;gto  "INC" 

180:  R[I+1]/V(2"N)*P;I[I+1]/V (2*N)*Z 

181:  R[J  +  1]*R[I+1]  ;I  [  J+l]  *I[I+1] 

162:  P*R[J+1] ;Z*I[J+1] 
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183: 

164: 

185: 

lb6: 

187: 

166: 

169: 

190: 

191: 

192: 

193: 

194: 

195: 

196  : 

197: 

198: 

199: 

200: 

201: 

2  02  : 

203: 

204: 

205: 

206  : 

207: 

206  : 

209: 

210: 

211: 

212 : 

213: 

214: 

215: 

216: 

217: 

218: 

219: 

220: 

221; 

222: 

222: 

224  : 

225; 

226 

227: 

228 

229; 

230; 

231 


"INC":R[I+1]  /y  (  2  Ai; ) -R  ( I +1]  ;  I  [1  +  1]  /v  (2*M*I  [1  +  1] 
"BB":next   I 

lr    rig  5; re i    11, I$,R[*] ;deg; ret 
a  eg  ;  r  e  c 

"SI"  :0+pi;pl-*p4 
ror    z=l    tc  p3 
p4/2-*p4 ;2*p2*p2 
ir    rrc(p4) #0;p2+l*p2 
int (p4) *p4 
next  Z ; ret 

"AIRY" : if    pl<0;beep;asp    "errcr-argument<0"  ; stp 
it    pl=0;  l-*r4;  ret    r4 
0->-r5;ii    pl>15;3mp    2 
2G+10*pl-pl~2/3-»-r6  ;  jir.p  -.2 
90  +  pl/2-*r6 

ir    pl<5;  6+pl-»-rl2;  jmp    2 
1 .4*pl+6  0/pl*rl2 
nr.ax  (int  irl2)  ,.int(3  +  pl/4)  )  +rl2 
tor    >:=rl2    tc    r6    fcy    3  ;le-2o  *r  6  ;  0*rl3-^r  14 
stg    10;  it   M/ 2=int(M/ 2) ;crg   lo 
fcr    J=l    tc  M-2;  2*  (K-J)-'*r8/pl-rl3+rl5;  r8*rl3 
rl5*r8;it    M-J-2=0;  r  15+r4 
'crr.f    10;  rl4+2*r8*£lgl0*rl4;  next   j 
2*rti/pl-rl3-rl5 
rl4+rl5+rl4;  r4/rl4+r4 

ir    acs  (r  4-r5).-at>£  (r  4*le-6)  <  =  0;  2*  r4/pl-r4  ;  r  et    r4 
r4*r5;next    h 

fceep;asp    "ACCURACY    NOT   OBTAINED" ; wait    1500 ;ret    r4 
"lA£LE":iat    pl,I$;R[*J 
tor    1=1    tc    32;rcr   J=l    tc    15 
wrt    "print. 2"  ,R[ 16 (1-1) +J] ; next  J 
wrt    "print.  3"  ,R[16I]  ;  next    I;  ret    G+Z 
"PLOT": la t    p2,I$,R[*j 
0*A;pl+B;min(R[  *]  )  *c  ;  max  (R  (  *  J  )  *D 
scl    A-.l  (B-A)  ,B-h.G5  IB-A)  ,C-.  l(D-C)  ,D-h.05  (D-C) 
a+E; 1G-*»F 
n    pl=512;64*G 
if    pl=256;32-G 
it    pl=64;6*G 
if    pl=32;4*G 
pit   B,C,1 

tor  l=£  tc  0  fcy  -G 
pit  I  (B-A)/E+A  ,C  ,2 
pit  I(B-A)/E+A  ,C+(u-C)/150,2 
pit  I(B-A)/E+A,C,2 
next  I 

tor  1=0  tc  F 
pit  A, I(D-C)/F+C,2 
pit  A+(6-A) /150,I(D-C)/F+C,2 
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2  22:  pit    A, I (D-C)/F+C,2 

22  3:  next    I; pen 

234:  csiz    1.2,1,.. 7 

235:  txa    1 

236:  tor    I=F    tc    0    ty  -1 

227:  pit   A-.075(B-A) ,1 (D-C)/F+C ,1 

238:  lbl    I/F 

229:  next    I 

240:  rxc   0 

241:  it    p3=0;"DATA    POINTS"-A$ ; 1-L 

242:  it    p3  =  2; "LINES/MICRCRADIAN"-A$ ; 1/(2*L) *L 

24  3:  tor    1=0    tc    E    by    G 

244:  pit    A+( I/E-.025J (B-A) ,C-.G25(D-C) , 1 

245:  lbl    I*L 

245:  next    I;  it    p3#G    ana    p3*2 ;  '\\ICRCRADIANS"*A$ 

24  7  :  it    p2  =  2;  S-*L 

248:  pit    .4 (B-A) +A A-.G5(D-C) +C ; lDl    A$ 

249:  pit  .-.07  (B-A)  +A  ,  ..3  (D-C  ) +C  ;  csiz    1.2,1, .7,90 

250:  "NORMALIZED    INTENSITY" +A$; lbl    A$ 

251:  csiz    ,5,l,1.5,0;0-»-I 

252:  pit    I  ,ir,ax  (R[*j  )  ,1 

253:  ror    1=1    tc   pl;pit    1-1 ,R[I ] ; next    I 

254:  csiz    1.2,1,,. 7,0 

255:  it    p2#0;gtc    +2 

256:  beep;ent    "PLCT    LABEL?" ,1$ 

257:  pit    .6B,.9D,1; lbl    "PLOT    CF" 

258:  pit    ,6B,.87D,l;lbl    1$ 

259:  pit    . 6B,.84D, l;lbl    ES,"    LASEF." 

26C:  pen;ctg     ;ret    0*Y 

*28629 
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0:  "COMPLIED    SOURCE    PATTERN": 

1:  ina    R;  10C0*A;r2*L-K:;  trk    l;fdf    10 

2:  for    1=1    tc    256 

3:  L*(L-])-F 

4:    F~2/(2*C~2)-G 

5:     if    G>13;jmp    2 

6:    A*exp(.-G)  *R[I]  ;I-*K;  jir.p   2 

7:    0-R[I] 

8:    next    I 

S:  "COMPUTED  SOURCE" -1$ 

10:  rcf  10,IS,R[*] 

11:    beep;ent    "PLOT   OF   COMPUTED    SOURCE  =  l   AND/OR    CCNT"  ,1 

12:    it    £«l;cll  "'PLOT' (32,10, Q) 

12:  "CONVERTS  9Si    TO  LSF": 

14:  Cll  'LSF" 

15:     "LSF    OF   COMPUTED    SOURCE" -1$ 

16:    rcf    10,1$, R[*j 

17:    ceep;ent    "PLOT    OF    LSF    OF    CCMP    SOURCE  =  l  AND/OR    CON.T"  ,Y 

18:    if    Y  =  l;cll  "'PLOT' (32,10,0) 

19:     "CALCULATION    OF    FXFORM   OF    LSF    OF    COMPUTED    SOURCE": 

20:    ina    I;ara    R*I 

21:    I[1]*R[1J 

22:    ror    1=2    tc    256 

22:    L[I] +R[I]+R[ 514-1 ]; next    I 

24:    R[256] -*R[257]  ;rcr    10,I$,R[*] 

25:    cll  "FXFORM ' 

26:     "FXFORM    OF    LSF    OF    C£LC    SOURCE" *I 5 ; ret    10,I$,R[*j 

27:     ceep;ent    "PLOT    OF    FXFM    OF    CALC    SRC£=1    AND/CR    CGNT" ,Y 

28:    it    Y  =  l;cll  "'PLOT" (256,10,2) 

29:     "CALCULATION    OF    PRODUCT   OF    2    FXFORMS": 

30:    ldf    10,I$,R[*] ;ina    I;ara    R+I;tr*    0 

31:     lat    10,I$,R[*] 

52:    tor    1=1    tc   256 

33:    R[I] *L[I]-R[I] ;next    I 

34:    rcf    10,I$,R[*] 

35:     "TRANSFER    FCN    OF    SOURCE*ATMCSPHERE*OPTTCS " : 

36:     lat    11 ,  1$  ,R[  *  J  ;  ina    I;ara    R-»I 

37:    1.(1]  *R[1]-R[1] 

38:    for    1=2    tc   256 

39:     R[I] *I[I]*R[I] +R[514-I] 

40:    next    I 

41:  R[256]-R[257] 

42:  "FXFM  OF  SOURCE*ATMOS*OPTICS" +1  $ 

43:  rcf  10,I$,R[*] 

44:  beep;ent  "PLOT  OF  FXFM  PRODUCTS  =  l  AND/OR  CON1"  ,Y 

45:  if  Y  =  l;cll  "PLOT"  (32,10,2) 
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"INVERSE    FXFORM    GIvES    TARGET    LSF": 

sfg    7;cll  "FXFCRM";cfg    7  ;  "INVERSE    FXFORM"*I$ 

rcf    10,I$,R[*] 

oeep;ent  "PLCT  OF  INV  FXFM=1  AND/OR  CCNT"  ,Y 

if  Y=l;cll  "'PLOT '(32,10,0) 

"CONVERTS  LSF  TO  PSF  BY  A EEL  TRANSFORM": 

Cll  "ABEL';  "ABEL  -XFORM  PSF"-I$ 

rci  10,1$ ,R[*] ; sfg  0 

beep;ent  "PLOT  OF  ABEL  XFCRM=1  AND/OR  COKT" ,Y 

if  Y=l;cll  "PLOT* '(32,10,0)  ;stp 

"CALCULATES  FRACTION  OF  POtvER  INSIDE  CIRCLE  OF  RADIUS  R" 

.25*7T*R[1]-R[1] 

for  1=2  tc  256 

2tt*I*R[I]+R[L-1]+R[I] 

next    I 

"POVvER    FRCTN    INSIDE    CIRCLE"  +1$ 

rcf    10,I$,R[*J 

ceep;ent    "PLOT    OF    POtoER=l  AND/CR    CONI"#Y 

if    Y=l;cil  "PLOT ' (32,1G,C) 

stp 

"LSF":ir    tlg5;lor    11  ,1  $  ,R  (  *]  ;nn.p    2 

ldf    10,I$,R[*] 

ina    I;ara   R>I;tor    1=1    tc    24 

l*J;csp    I;I[I]+Q 

V  (I*I+J*J)+R 

2*  ( (l-frc(R)  )*I[int(R)]+frc(R)*L[int{R)+l]  )+Q*g 

J+l+J;  if    R<24  ;  jlT.p    -2 

y*i[ij 

next    I; ina    R;ara    I+R; if    rlgb;rcr    11,1  $ ,R[*]  ;ctg    ;ret    pi 
ret   pi 

"FXFORM"  : rad ; 9+N ;  ina    I;  if    flg5;laf    11 , 1  $  ,R  [*]  ;ijmp    2 
lot     1U,I$,R[*] 
tt/2"(N-1)  +T 

for    M=l    tc  N;  2~ (K-M) *rG 
for    J  =  0    tc    2"(ML-l).-l;cll  "Bl"(J  ,P  #N-1) 
.ccs(F*T)*C;sin(P*T) * (1-2* tig  7 ) -P 
for    I=2*rO*J+l    tc    2*rQ*J+rO 
R[I]*rl;R[I+r0] -r2 
I [IJ  +  r3; L[I  +  rCJ  -r4 

rl+r2*C+r4*P*R[I] ; r 3+r4*C- r2*P+I [I] 
rl-r2*C-r4*P-*R[  I+r0]  ; r3-r4*C+r 2*P+I [I+r0 ] 
next    I; next   J ;asp    K; next  m 
for    1  =  0    tc    2"N-l;cll  "BI"(I  ,J,N) 
if    I-J>0;gto    "BB" 
if    I=J;gto    "INC" 

R[I  +  l]/V(2'%N)*P;L[I+l]/v(2*ti)+Z 
R[J+1]+R[I+1]  ;I  [J+l]  *I[I  +  1] 
P*R[J+1] ;Z*I [J+l] 
,,INC":R[I+1]/V(2"'N)*R[I+1]  ;  I  [I  +  l]/v  (  2*N)  +1  [  1  +  1] 
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95:    "EB" :next    I 

96:    cegjret 

97:     "BI":0+p2; p]*p4 

98:    for    Z=l    tc   p3 

99:    p4/2+p4;2*p2+p2 

100:    if    rrc(p4) #0;p2+l+p2 

101:    int(p4)-p4 

102:    next   Z;ret 

103:     "A£EL":ina    R;lcf    10  ,1  $  ,R  [*  ]  ;  R  [  1]  -N 

104:    1.4*R[l]-1.8*R[2]-h.4*R[3]*R[l] 

105:     for    L=2    tc    64 

106:    R[I]-M 

107:     .4*N+.2*M-.6*R[I+1] -R  [  I  ] ;M-N;next    I 

108:    for    1=1    tc    64 

109:    R[I]/(2*V  (  (I-h.l)  "2-1*1.)  )  +R  [  I  ] 

110:    for    J=I+1    tc    64 

111:    R[I]+R(J]/V  (  (J-h.l)  "2-I*I)*R[I] 

112:    next  J 

J13:    if    R[L]<.01;G+R[I] 

114:    R[I  ]/tt*R[L]  ;asp    I;  next    I 

115:    for    1=65    tc    512; G+R[I] ; next    I; ret 

116:     "PLOT" :ld£    p2,I$,R[*] 

117:    0-"A;pl-B;min(R[*]  )  *C;  max  (R[*j  )  +C 

116:    scl   A-.l(B-A)  ,B-h.05  IB-A)  ,C-.  l(D-C)  ,D+.05(D~C) 

119:    B+E;10+F;if    f IgO ; pr t   A  ,B ,C,D 

120:     if    pl=512;  64-*G 

121:     if    pl=256; 32-G 

122:     if    pl=64; b*G 

123:    if    pl=32;4*G 

124:    pit   5,C,1 

125:    for    l=E    tc    0    by    -G 

126:    pit    I (E-A)/E+A  ,C  ,  2 

127:    pit    I(E-A)/E+AfC+(D.-C)/150,2 

128:    pit    I (B-A)/E+A,C,2 

129:    next    I 

130:    for    1=0    to   F 

131:    pit   A,I(D-C)/F+C,2 

132:    pit    A+(B-A)/150,I (D-C)/F+C,2 

133:    pit    A,I(D-C)/F+C,2 

134:    next    I; pen 

135:    csiz    1 . 2, 1,  .7 

136:    fxa    1 

137:    for    I=F    to    0    by   -1 

138:     pit    A-.075(B-A) ,1 (D-C)/F+C,l 

139:     lbl    I/F 

140:    next    I 

141:    fxd    0 

142:     if    p3=2; HLINES/MICRCRADIANM*A$ ; 1/(2*L) +L 

143:    for    1=0    tc   E    by   G 
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144:  pit    A-r  (  i/E-  .  025.)  (B-A)  ,C-.025{D-C)  ,1 

145:  lbl    I*L 

l4o:  next    Ijif    p3#2;  ,,MICRORADIANS"+A$ 

147:  it    p3=2;S-L 

148:  pit    .4(B-A) +A,.-.05(D-C)+C;  lbl   A$ 

149:  pit  .-.07  (B-A)  +A  ,  ..3  (D-C) +C  ;  csiz    1.2,1,  .7,90 

150:  "NORMALIZED    INTENSITY" +A$; lbl    A$ 

151:  csiz    .5,1,1.5, G;0+I 

152:  pit   I ,max(R[*J ) ,1 

153:  for    1=1    to   pi; pit    I-1,R [ I ] ; nex t    I 

154:  csiz    1.2 ,1*. 7,0 

155:  pit    .6£,.SD,l;lol    "PLOT   OF" 

156:  pit    .6B,  ,87D,1; lol    1$ 

157:  pit    .6B,  .84D, l;lbl    E$,"    LASER" 

158:  "REI" :pen;ctg     ;ret    0-i 

*  1 7  5  2h 
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